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The analysis of the resonant spin-flavour (RSF) solutions to the solar neutrino problem in the framework of 
simplest analytic solutions to the solar magneto-hydrodynamics (MHD) equations is presented. We performed the 
global fit of the recent solar neutrino data, including event rates as well as day and night recoil electron spectra 
induced by solar neutrino interactions in SuperKamiokande. We compare quantitatively our simplest MHD-RSF 
fit with vacuum oscillation (VAC) and MSW-type (SMA, LMA and LOW) solutions to the solar neutrino problem 
using a common well-calibrated theoretical calculation and fit procedure and find MHD-RSF fit to be somewhat 
better than those obtained for the favored neutrino oscillation solutions. We made the predictions for future 
experiments to disentangle the MHD-RSF scenario from other scenarios. 



1. Introduction 

The problem of disagreement between solar 
neutrino data and theoretical expectations has 
been a long-standing problem in physics. The 
most popular solutions of the solar neutrino 
anomalies have as a basis the idea of neutrino os- 
cillations, either in vacuum or in the Sun due to 
the enhancement arising from matter effects [Q. 
There are alternative interpretations the problem. 
Here we will re-analyze the status of resonant 
spin-flavour solutions to the solar neutrino 
problem in the light of the most recent global set 
of solar neutrino data. In contrast to previous 
attempts ||-|^] we will adopt the general frame- 
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work of self-consistent magneto-hydrodynamic 
(MHD) models of the Sun || . For definiteness we 
will concentrate in the recent proposal of Ref. || . 
We perform global fits of solar neutrino data for 
realistic solutions to the magneto-hydrodynamics 
equations inside the Sun. This way and by ne- 
glecting neutrino mixing we obtain the simplest 
MHD-RSF solution to the solar neutrino prob- 
lem, characterized by two effective parameters, 
Am 2 and i*ivB± max , B± max being the maximum 
magnitude of the magnetic field inside the con- 
vective region. We find that our simplest two- 
parameter MHD-RSF fits to the solar neutrino 
data are slightly better than those for the oscil- 
lation solutions. The required best fit points cor- 
respond to maximum magnetic field magnitudes 
in the convective zone smaller than 100 KG. We 
briefly discuss the prospects to distinguish our 
simplest MHD-RSF scenario from the neutrino 
oscillation solutions to the solar neutrino prob- 
lem. 
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2. Static Magnetic Field Profiles in the 
Sun 

In solar magneto-hydrodynamics Q (MHD, for 
short) the corresponding magnetic field profiles 
are rather complicated and difficult to extract. 
However, there are stationary solutions which are 
known analytically in terms of relatively simple 
functions ||. 

We consider the magnetic field profile which 
are only solutions to the equation for a static 
MHD plasma configuration in a gravitational 
field, given by 
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Vp j x B + pV$ = 0. 
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This static MHD equations correspond to a quiet 
Sun and they admit axially symmetric solutions 
in the spherically symmetric gravitational field 
which can be simply expressed in terms of spher- 
ical Bessel functions and were first discussed in 
Ref. |J. The model magnetic field depends on z k , 
the roots of the spherical Bessel function / 5 / 2 = 
\fzJ b / 2 { z )- Taking into account the boundary 
condition that B vanishes on the solar surface the 
magnetic field have the analytical form 
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Figure 1. The perpendicular component of B for 
various fc-values 1 (solid), 3 (dashed) and 10 (dot- 
ted). 



The distance r has been normalized to Rq = 1. 
In our calculations we have averaged over polar 
angle 9. 

Parameter B core is a central magnetic field. In 
Fig. 1 we display the perpendicular component 
of B for various fc-values 1, 3 and 10, which cor- 
respond to the roots z\ = 5.7, 23 = 12.3 and 
zio = 34.5, respectively. 

We have to discuss the astrophysical restric- 
tions on the free parameters B core and k charac- 
terizing the model. The magnitude of a magnetic 
field at the center of the Sun is constrained by the 
Fermi- Chandrasekhar limit Q which requires an 
upper bound on B core < 2 MGauss. 

The possible values of k can be constrained 
by taking into account that in order to justify 
the use of a stationary solution, it is necessary 
that the diffusion time due to ohmic dissipation 
must be less than solar life time. The simple 
estimations give us that reasonable values of k 
are less 10. It is commonly accepted that mag- 
netic fields measured at the surface of the Sun 
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Figure 2. Magnetic field configurations obtained 
by combining individual modes for different Um 
values, 5, 6 and 10 (dotted, dashed and solid). 
Summing up to higher modes achieves better lo- 
calization of the field in the convective region 
(solid). 



are weaker than within the convective zone in- 
terior where this field is supposed to be gener- 
ated. On the other hand the general knowledge 
of the solar magnetic field models is that the 
magnetic field increases at the overshoot layer, 
while being small at the solar interior, a picture 
rather opposite to the one we have seen in Fig. 1. 
The correct way is to use the linear nature of 
the basic equilibrium MHD equation in eq. (Q). 
This implies that any linear combination of solu- 
tions B k (k — 1,2, ... , fcjw, for some fixed number 

k M < 10) B = c\B\ + c 2 B 2 + ... + c M B M is also 
a solution. We will require that combined mag- 
netic field is equal to zero in the center of the 
Sun and it's total energy must be minimal in the 
region below the bottom of the convective zone, 
characterized by a certain value of Tq. 

The procedure sketched above provides a con- 
sistent method for combining individual mode so- 
lutions Bk of the static MHD equation (Fig. 2). 



3. Fitting the Solar Neutrino Data 

We will neglect neutrino mixing and consider 
the case of active-active neutrino conversions. In 
this case the v e — > vg_ conversions are described 
by the master Schrodinger evolution equation 



V e -5 



p v B^ 



, (4) 



where /i„ denotes the neutrino transition mag- 
netic moment || in units of 10~ n /is, £ denot- 
ing either /i or r. Here B± = B x ± iB y and 
5 = Am 2 /4-E is the neutrino mass parameter; 
V e (t) = G F V2(p(t)/m p )(Y e - Y n /2) and V e {t) = 
GFV2(p(t)/m p )(—Y n /2) are the neutrino vector 
potentials for v e and Vi in the Sun given by the 
abundances of the electron (Y e — m p N e (t) / p(t)) 
and neutron (Y n — m p N n (t) / pit)) components. 
In our numerical study of solar neutrino data we 
adopt the Standard Solar Model density profile of 
ref.flO). 

We solve Eq. (||) numerically by finding a solu- 
tion of the Cauchy problem in the form of a set of 
wave functions v a (t) = \ v a (t) \ e 1 *" ^ from which 
the neutrino survival probabilities P aa {t) — v*v a 
are calculated. They obey the unitarity condition 
J2a = 1 where the subscript a denotes a = e 
for v e and a = i for vi respectively. 

The typical neutrino survival probability P ee 
calculated in the MHD-RSF scheme from eq. (|4|) 
plotted versus E/Am 2 in Fig. 3. 

We determine the allowed range of oscilla- 
tion parameters using the total event rates of 
the Chlorine [O, Gallium and Super- 

Kamiokande ]15| , |16[ (corresponding to the 1117 
days data sample) experiments. For the Gallium 
experiments we have used the weighted average of 
the results from GALLEX+GNO and SAGE de- 
tectors (see Table [[]). We have also included the 
Super-Kamiokande electron recoil energy spec- 
trum measured separately during the day and 
night periods. For details on the statistical anal- 
ysis applied to the different observable we refer to 
Ref. . 

We have found that allowed regions of neutrino 
parameters are pretty stable and does not depend 
significantly on the choice of kM and ro allowed 
by astrophysics. In Fig. 4 we display the region of 
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Table 1 

Solar neutrino rates measured in the Chlorine, Gallium and Super-Kamiokande experiments. 



Experiment 



Rate 



Ref. 



Units 



R 



BF2D0D 



Homcstakc 2.56 ±0.23 

GALLEX+GNO+SAGE 74.66 ± 5.2 

Super-Kamiokande 2.40 ± 0.08 



12 -$A 
16] 



SNU 
SNU 

10 6 cm~2 s" 



7.8 ± 1.1 
130 ±7 
5.2 ±0.9 




E/Am 2 (eV~ 1 ) 



Figure 3. Typical MHD-RSF neutrino survival 
probability P ee versus E/Am 2 . 



MHD-RSF parameters allowed by the solar neu- 
trino rates for the case Um = 6 and r = O.6i?0. 
We can see that there are several allowed regions 
for different values of the magnetic field. 

Apart from total event rates the water 
Cerenkov experiment also measures the zenith an- 
gle distribution of solar neutrino events as well as 
their electron recoil energy spectrum with their 
recent 1117-day data sample Jl6) . The predicted 
spectrum is essentially flat except for the upper 
part of the Am 2 region. As an example, we show 
in Fig. 5 the excluded region at 99 % CL for the 
case fcjvf = 6 and ro = 0.6. For this reason, the 
allowed regions are slightly modified by the inclu- 
sion of the energy spectrum data for the day and 
night periods. In Fig. 6 we have presented the 
results of global fit analysis. 



£~J0" 7 I . 
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B Mfp (kG) 



Figure 4. MHD-RSF 90% CL (light) and 99% 
CL (dark) regions of Am 2 versus Bi max (KG) al- 
lowed by the rates given in table |l|, for r = 0.6 
and k,M = 6. 



In the case of active-sterile MHD-RSF conver- 
sions we obtain that, the rates fit Xrates ^ s worse 
than for the active-active case. 

4. Discussion &: Conclusions 

From the results of the previous section it fol- 
lows that our MHD-RSF solution to the solar 
neutrino problem provides a good description of 
the most recent solar neutrino data, including 
event rates as well as zenith angle distributions 
and recoil electron spectra induced by solar neu- 
trino interactions in Superkamiokande. We have 
shown that our procedure is quite robust in the 
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Table 2 

Best fit points and the corresponding probabilities for different solutions to the solar neutrino problem jjj 
The top row corresponds to the MHD-RSF solution presented here. 



Solution 


Am 2 


B _Lmax 


Xtnin %) 




MHD - RSF a 


1.1 x 1CT 8 


80 


30.4 (73) 


This work 


MHD - RSF S 


1.1 x 1(T 8 


77 


34.9 (52) 






Am 2 

c n „ m-6 


sin 2 (2??) 

no., lr>-3 


X 2 mm (Prob %) 


Ref. 



LMA 
LOW 
QVO 
VAC 
SMA, 



3.2 x 10~ 5 
1.0 x 10~ 7 

2.3 x 10~ 9 
6.7 x 10~ 10 
3.9 x 10~ 6 



0.75 
0.93 
0.96(d) 
0.93(d) 
2.4 x 10- 



33.4 (59) 
37.4 (40) 
40.3 (29) 
39.3 (32) 
39.6 (31) 



[LZI 

[17| 
[17| 
[17| 



no-osc 



91 (3 x 10~ 4 ) 



sense that the magnetic field profile has been de- 
termined in an essentially unique way. This ef- 
fectively substitutes the neutrino mixing which 
characterizes the oscillation solutions by a single 
parameter B± max characterizing the maximum 
magnitude of the magnetic field inside the con- 
vective region. The value of Hm characterizing 
the maximum number of individual modes su- 
perimposed in order to obtain a realistic profile 
and the parameter r$ characterizing the location 
of the convective region are severely restricted. 
The allowed fcjvf values are restricted by ohmic 
dissipation arguments to be lower than 10 or so, 
while ro is close to O.6i? . We have found that 
our solar neutrino fits are pretty stable as long 
as fcjw exceeds 5 and ro lies in the relevant nar- 
row range. Therefore our fits are effectively two- 
parameter fits (Am 2 and B± max ) whose quality 
can be meaningfully compared with that of the 
fits obtained for the favored neutrino oscillation 
solutions to the solar neutrino problem. 

In table || we compare the various solutions of 
the solar neutrino problem with the MHD-RSF 
solutions for the lower magnetic field presented 
here. 

Clearly the MHD-RSF fits seem somewhat bet- 
ter (though not in a statistically significant way) 
than those obtained for the MSW effect jl7| as 
well as just-so solutions |]| . 

We determined the expected solar neutrino 
rates at SNO within the framework of our MHD- 



RSF solution to the solar neutrino problem. We 
used the cross sections of the CC and NC vd reac- 



tions given by ref. 19 and the best-fit points we 
have determined in the present paper. For def- 
initeness we have considered the global best fit 
points and local minima for B± max < 100 KG for 
the case A:m = 6 and ro = 0.6 and active-active 
MHD-RSF conversions. 

We have calculated the neutral-to-charged- 
current event ratio (NC/CC for short) and our 
results are presented in Fig. 7. 

Clearly from Fig. 7 we see that there is a sub- 
stantial overlap between our MHD-RSF predic- 
tions and those found for each of the oscillation 
solutions (SMA, LMA, LOW, VAC). The over- 
lap is especially large between the LMA and the 
MHD-RSF solutions. Taking into account the 
present theoretical uncertainties and a reasonable 
estimate of the experimental errors attainable, it 
follows that an unambiguous discrimination be- 
tween our MHD-RSF solution and the neutrino 
oscillation-type solutions to the solar neutrino 
problem on the basis of the averaged event rates 
seems rather difficult. 
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Figure 5. MHD-RSF 99% CL regions of Am 2 ver- 
sus B± max forbidden by the recoil electron spec- 
trum data for ro = 0.6 and Hm = 6. 



Figure 6. MHD-RSF 90% CL (light) and 99% 
CL (dark) regions of Am 2 versus B± max (KG) al- 
lowed by the global solar neutrino fit, for ro = 0.6 
and /cm = 6. 
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